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Previous studies jlj have provided strong support for a local, oscillatory approach to the singu- 
larity in U(l) symmetric, spatially inhomogeneous vacuum cosmologies on T 3 x R. The description 
of a vacuum Bianchi type IX, spatially homogeneous Mixmaster cosmology (on S 3 x R) in terms of 
the variables used to describe the (7(1) symmetric cosmologies indicates that the oscillations in the 
latter are in fact those of local Mixmaster dynamics. One of the variables of the U(l) symmetric 
models increases only at the end of a Mixmaster era. Such an increase therefore yields a qualitative 
signature for local Mixmaster dynamics in spatially inhomogeneous cosmologies. 
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I. INTRODUCTION 

o , 

In the 1960's, Belinskii, Khalatnikov, and Lifshitz (BKL) developed a formalism to analyze the approach to the 
I | singularity in spatially inhomogneous cosmologies j^-Qj- They argued that, while one can construct, over any short 
time interval, a local Kasner solution [Q , this type of solution cannot be maintained indefinitely since terms inconsistent 
H ' with it cannot remain subdominant. Rather, the generic solution will behave as a spatially homogeneous Mixmaster 
, cosmology Q at every spatial point. For a long time, this analysis remained controversial, primarily due to its 
■ non-rigorous nature. 

Recently, in a series of papers [[I]J^-|l|, we and our collaborators have used an analytic method, simpler than but not 
unrelated to that of BKL, along with detailed numerical simulations to provide strong support for the BKL picture. 
In addition, our methods have clarified how Einstein's equations control the approach to the singularity. The essence 
of our method of consistent potentials (MCP) developed originally by Grubisic and Moncrief [Q is to solve the 
velocity term dominated (VTD) equations obtained by neglecting all terms containing spatial derivatives in Einstein's 
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equations. We then substitute the asymptotic form of the VTD solution (with spatially dependent parameters) into 
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X 



the neglected terms. If these terms are all exponentially small, the VTD solution is consistent. This in turn supports 
a conjecture that the approach to the singularity is asymptotically VTD (AVTD) — i.e., at almost every spatial 
point, the full solution comes arbitrarily close to a VTD solution as the singularity is approached On the other 
hand, if one or more of the neglected terms grows exponentially, then the VTD solution is inconsistent so that the 
conjectured behavior is not AVTD. If two or more previously neglected terms will alternately always grow for any 
range of VTD solution parameters, then the MCP predicts that the solution will be oscillatory and the oscillations 
will persist forever. If there exists a range of parameter values consistent with a VTD solution almost everywhere, 
the prediction is that oscillations will persist only until the parameters are driven into the consistent range lfjl. For 
the MCP to predict the correct behavior, the only dynamically significant time dependences, sufficiently close to 
the singularity, must be the linear ones from the VTD solutions which appear in the arguments of the exponential 
potentials. Perhaps surprisingly then, the MCP predictions have been verified by numerical simulation of the full 
Einstein equations in all models we have studied so far fO] . 

The previous discussion of the application of the MCP makes no appeal to an invariant description of the system. 
One then wonders whether or how this picture depends on the choice of variables and of spacetime slicing. There 
is also an issue of the nature of oscillatory behavior observed in spatially inhomogeneous cosmologies |Q. The BKL 
claim is that these models should evolve as a vacuum Bianchi type IX (Mixmaster) universe at every spatial point. In 
the MCP applied to the minisuperspace (MSS) picture of vacuum, homogeneous Bianchi type IX models, the VTD is 
the Kasner solution. This solution is inconsistent with maintaining exponential smallness of the three exponential wall 
terms in the MSS potential — leading to the usual infinite sequence of bounces off a closed triangular potential || . 
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However, we found that vacuum U(l) symmetric cosmologies on T 3 x R exhibit oscillations between two potentials 
(rather than three) at every spatial point At the time, we did not know if these represented "true" Mixmaster 
oscillations or something else. 

In this paper, we shall describe spatially homogeneous Bianchi type IX universes on S 3, x R as U(\) symmetric 
cosmologies with the same topology by rewriting the Bianchi type IX spatial metric in a coordinate frame Jl6| . The 
objective of this exercise is to see how the standard Mixmaster oscillations are described in the U(l) variables used in 
our previous studies. We find that the relationship between the MSS and ?7(1) descriptions of the Mixmaster universe 
is not trivial. In the coordinate representation of the Bianchi type IX homogeneous space, one of the Killing vectors 
is still manifest and is taken to play the role of the spatial U(l) symmetry. The norm of that Killing vector is the 
key variable in our previous description of generic, spatially inhomogeneous U(l) models. Here we find its behavior 
to be dominated by the largest logarithmic scale factor (LSF) of the standard BKL description. The bounces off one 
of the U(l) potentials occur when the time derivative of the dominant LSF changes sign. This is what happens in the 
standard MSS picture in a bounce off any of the three potential walls. The other U(l) potential causes a bounce when 
one LSF loses its dominant place to another in the MSS picture. This correspondance will be discussed extensively 
in this paper. 

Although numerical analysis of inhomogeneous U(l) symmetric cosmologies revealed only the two types of bounces 
described above, our study of the {/(l)-MSS correspondance for the Mixmaster universe demonstrates that a third 
type of bounce can occur — at the end of a BKL era jj| — when the monotonically decreasing smallest LSF starts 
increasing and overtakes the middle one. This yields a qualitative prediction — not yet observed numerically — which 
can distinguish between local Mixmaster dynamics and other oscillatory behavior. 

In Section II, we shall review the MSS picture and, in Section III, the U(l) variables. Note that in this paper, we 
shall often use "{/(l)" to refer to properties of the U(l) symmetric cosmologies. The correspondance between the 
U(l) models and the Mixmaster models will be discussed in Section IV with the implications of the correspondance 
given in Section V. Conclusions will be drawn in Section VI. 

II. THE MSS PICTURE OF MIXMASTER DYNAMICS 

A vacuum, diagonal Bianchi type IX Mixmaster cosmology is described by the metric jl7| 

ds 2 = ~A 2 B 2 C 2 dr 2 + A 2 {a 1 ) 2 + B 2 (a 2 ) 2 + C 2 (a 3 ) 2 (1) 
where the scale factors A, B, and C are functions of r only and 

da*= l -e) k aifxa k . (2) 

It is often convenient to define the LSF's a, £, and 7 by 

A = e a , B = e c , C = e 7 . (3) 

The SU(2) symmetry (§ may be realized in a coordinate frame by 

a 1 = cos 4>d9 + sin 8 sin (p dtp, 
a 2 = — sin 4>d9 + cos <j> sin 6 dip, 

a 3 =d(p + cos 6 dip (4) 



on S . The choice of time coordinate 



N dr = ABC dr = dt, (5) 



where t is comoving proper time, is that used by BKL. 

An alternate, but equivalent description, was given by Misner || in terms of the MSS variables O, the logarithmic 
volume, and p)±, the anisotropic shears, where 

a = n - 2/3+, 

C = « + /3+ + V3/3_, 

1 = n + [3 + - V3/3_. (6) 



2 



This variable choice conveniently allows Einstein's equations to be obtained by variation of the superhamiltonian 
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(7) 



where g is the determinant of the spatial metric gtj with conjugate momenta 7r y while 3 R is the scalar curvature of 
gij and H = is the Hamiltonian constraint. With Pq, p± canonically conjugate to CI, (3± respectively, and the time 
coordinate choice N = ^/g/sinfl, Eq. (Q) becomes 



for the MSS potential 



Vix 



2H = = -p 2 n +p\+ p 2 _ + Vi X (p±,Sl) 



(8) 
(9) 



When Vix is exponentially small, the metric is locally (neglecting the topology) the Kasner solution. Generically, in 
the evolution toward the singularity, one of the first three terms on the right hand side of Eq. (||) — the one associated 
with the largest LSF — will grow. A "bounce" off the exponential potential will change the sign of the time derivative 
of the dominant LSF changing one Kasner solution into another. Conservation of momentum in the bounce was first 
used by BKL to relate the two asymptotic Kasner solutions. In terms of the LSF's, Einstein's evolution equations for 
the Mixmaster model are (for overdot indicating d/dr) 



a 
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If the initial Kasner solution is characterized by a given {a, C>j} an d the dominant potential term is e 4a , then from 
Eqs. ( |l0| ) the final asymptotic Kasner solution with {a', 7'} will have a' = —a, (' — C + 2d, 7' = 7 + 2a. Various 
schemes exist to encode this information [B 18 13]. If, initially, a > ( > 7, then initially a > while ( < 0, 7 < since 
the (collapsing) Kasner solution generically has only one expanding direction. This means that a increases while the 
others decrease. After a typical bounce, the above bounce rules yield a' < 0, £' > while 7 has increased but is still 
negative. In subsequent bounces, a and £ will continue to oscillate while 7 decreases monotonically but with ever 
decreasing |7|. (Each Kasner regime between the bounces is called an epoch.) Eventually, 7 will change sign. This 
is called the end of an era. If a > £ > 7 prior to this era ending bounce, we find that, asymptotically, 7' > £' > d' 
after the bounce. Subsequently, ^ and 7 oscillate while a decreases monotonically. This sequence of eras apparently 
pQ-E2[ persists indefinitely. Part of a typical Mixmaster evolution is shown in Fig. 1. 



III. U(l) SYMMETRIC COSMOLOGIES ON S :i x R 



U(l) symmetric cosmologies on S 3 x R are described by the metric 



ds 2 



~ 2 ^{-N 2 dr 2 + g ab (dx a + N a d T ){dx b + N b dr)} + e 2v (d^ + cos8d(j) + f3 a dx a + drf 



(11) 



where the symmetry direction is tp, and the other spatial directions are {x a } — {9, </>}. The metric variables are 
assumed to be functions of 9, <fi, and r. The norm of the Killing field is e v , (3 a are the "twists", e 2A is the determinant 
of the 2-metric g a b = e A e a b and e a b is parametrized by x and z via 
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Note that the metric ( |Tl| ) differs from that for T 3 spatial topology given in |lj . It is convenient to make a canonical 
transformation from the twists and their conjugate momenta e a to the twist potential uj and its conjugate momentum 
r. It is also convenient to define the spacetime slicing by zero shift and lapse N sin 9 = \fg = e A where g is the 
determinant of the 2-metric g a b. Einstein's equations are found by variation of 
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+ { (e V 6 ) , ab ~ (e A e ab ) , a A, 6 +e A [(e~ 2z ) ,„ a:, - (e" 2 *) ,„ ] 

(13) 



+2e W,a <£,6 +ie A e- 4 ^e Q& W , a u, b 

where Ti. = is the Hamiltonian constraint and the overall trigonometric factor comes from N / yfg. 

An excellent approximation to the behavior seen in numerical simulations of these models (albeit on T 3 rather than 
S 3 ) is that at each spatial point a Kasner-like phase characterized by Q| 

If = —v v t, p = — 4u,p, 2 = —v z t, 

Pz = -4t!a, A = -W A T, p A = "A, 

w = cj , r = r , a; = x , Px=P° (14) 

(where t> v , t> 2 , «a, wqj foj ^Oj an d are functions of 8 and <j> but independent of r) is followed by a "bounce" off one 
of the three potentials 



Vi = 
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2 P«e ■ 



(15) 

In Eqs. (|l5|), the potentials are defined to be identical to those in JjJ. (It is possible to incorporate the prefactor N/y/g 
from Eq. (|13[) into the potentials. However, this factor depends on topology. In Section V we shall see that only the 
(subdominant) logarithm of this factor would appear in the variables of interest and would not affect the qualitative 
behavior.) If one assumes incident and final asymptotic solutions of the form ( |l4| ) at a given spatial point and that r, 
lu, x, and p x are constant in time through the bounce, the remaining momenta p, p z , and pa at that point will change 
during a bounce off one of the potentials according to the following rules: 

(1) For a bounce off V±, p — > — p with the other momenta conserved. 

(2) For a bounce off V%, 

P~> -P-Pz- 2pA, 

Pz -> - (~2p + p z - 2p A ), 

PA^^(2p + Pz + 6 PA ). (16) 

These rules are obtained by recognizing that, for z large and negative, e ab is dominated by e~ 2z so that V2 ~ 
j( e A-2z-4ip w j iere k i s approximately constant in time. Then A — 2z — Aip defines the direction (in a local MSS) in 
which the conjugate momentum changes sign during a bounce with momenta in the plane orthogonal to that direction 
conserved. 

(3) For a bounce off V3, p z — > — Pz with the other momenta conserved. 

(4) In all cases, the momenta must satisfy an asymptotic form of the Hamiltonian constraint ( fL3"| ) given by 

-Hum = -\pl + \p 2 + \p 2 z=V- (17) 

Eq. ( JjL7\) is obtained from Eq. (|l^) by assuming that all terms containing exponential potentials are exponentially 
small. 

In numerical simulations performed so far |jj , the dominant local behavior is oscillation of ip due to bounces off Vi 
and Vi. A bounce off V3 has been seen only if p z > initially. Note that the rule ( |l6| ) for a bounce off V2 becomes 
p — > —p if \p/p z \ << 1 implying p z 1=3 2pA (from Eq. (|l7|)). Our simulations provide strong support for BKL's claim 
that the approach to the singularity in generic cosmological spacetimes is oscillatory. The remaining question then is 
whether or not this oscillatory behavior represents local Mixmaster dynamics. 
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IV. THE MIXMASTER UNIVERSE AS A (7(1) SYMMETRIC COSMOLOGY 



By writing the Mixmaster metric (|l[) (denoted g IX in this section) in the coordinate frame on S 3 given by Eq. (Q) 
and comparing it to the U(l) metric ( |ll|) (denoted g u in this section), each U{1) variable ip, z, A, u>, or x and its 
conjugate momentum may be expressed in terms of the BKL scale factors A, B, and C and their time derivatives. 
The explicit Killing direction ip will be identified with the U(l) symmetry direction x 3 so that e 2ip = or 

e 2ifi = A 2 sin 2 9 sin 2 + B 2 sin 2 6 cos 2 + C 2 cos 2 0. (18) 

The twists /3# and in g^ 7 are determined from e 2v (3g — g 1 ^ and e 2lp (f3<p + cos 9) — g 1 ^ so that 

(A 2 — B 2 ) cos 4> sin 6* sin 



fio 



C 2 cos 2 6» + sin 2 9(B 2 cos 2 + A 2 sin 2 . 
C 2 cos<9 



9 C 2 cos 2 6» + sin 2 9(B 2 cos 2 + A 2 sin 2 0) 
Then the U{\) model 2-metric (here denoted as g a b = e A e a fc) is found from 



cos6>. (19) 



e 



- 2 *g ee + e 2 ^ 2 , 



9H = e ~ 2ip 9e<p + e 2ip (3 e {^ + cos 9), 

g'4 = e- 2 ^ + e 2 ^(/3 + cos6') 2 . (20) 

Comparison of Eqs. (|l|) and (|ll|) with Eqs. (|l^) and (|l^), computation of the determinant of g, and the expression of 
e a b in terms of x and z (see Eq. (12h) leads to 



ABC sin9 J C 2 cos 2 6> + sin 2 9(B 2 cos 2 </> + A 2 sin 2 0), (21) 



AABC sinO J C 2 cos 2 6> + sin 2 9(B 2 cos 2 + A 2 sin 2 0) 

P 2 ^ = V (22) 

A 2 B 2 +A 2 C 2 +B 2 C 2 ~A 2 B 2 cos(29) + (A 2 -B 2 )C 2 cos[2(9-(j))Y V ' 



A 2 B 2 + [5 2 C 2 + A 2 (C 2 - B 2 )] cos(26>) + {A 2 - g 2 )C 2 cos(20) 
A 2 S 2 + A 2 C 2 + B 2 C* 2 - A 2 B 2 cos(26») + (A 2 - B 2 )C 2 cos[2(<9 - 0)] 



(23) 



We note that 



F g = e A = ABCsin9e v . (24) 



Our gauge condition in jy (modified for S* 3 ) is = e A / sin^. But this is the "2-dimensional lapse." The 3-dimensional 
lapse is N = Ne~ v so that 

N = ABC (25) 

which is just the usual BKL lapse condition ||. This means that the time variable r chosen in £7(1) models is the 
same as the usual BKL time r. 

It is now necessary to consider the canonical transformation which replaces /3 a and its conjugate momentum e a 
subject to the constraint |f23[ 

e a , a -0 (26) 
by to and its canonically conjugate momentum r where 

e = e UJ, b (27) 
for e ab — —e ba . In these variables, e° may be found from the Einstein equation for the time derivative of j3 a |23j|: 
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(28) 



Wc find 



( = S \ r — {BCA sin <fi cos <j) — ACB sin <fi cos < 



=4> 



N 

2 sin 9 cos 9 



u,4> ■ 



N 



- {ABC - BCA sin 2 - ACB cos 2 



(29) 



to give 



sin 2 9 



-{-ABC + BCA sin 2 (j) + ACB cos 2 <j>) + k{r) 



where k(r) is an as yet undertermined function of time. 

The momentum conjugate to cj is defined by an integrability condition [^3| so that 



(30) 



r = f3e,<p-f3<p,8 + sin9 
= e-^ {{A 2 - B 2 ) sin 3 9{B 2 cos 2 <j> - A 2 sin 2 0) 

+ C 2 sin 6»[A 2 + B 2 -C 2 + sin 2 9{C 2 - A 2 cos 2 <f) — B 2 sin 2 <f>)]} 



(31) 



It is easily checked that J s2 r — An as is required. From the U{1) symmetric model Hamiltonian (13), we obtain 



uj, t = re Av 



N 

71 



(32) 



so that 



fc(r) = - 



A B 
A + B 



ko 



(33) 



where fco is a constant. 

The other momenta — p, p z , p\, p x — may also be found from the equations of motion — for ip, T) z, T , A, r , x, T 
respectively and are given for completeness in the Appendix. 

Given all the variables, it becomes a simple matter to construct all the terms in the U(l) Hamiltonian (|l3|). In 
particular, we shall consider 



Vi 



1 



2 4uj 

r e v 



{ {A 2 - B 2 ) sin 3 9 (B 2 cos 2 - A 2 sin 2 cf) 

+C 2 sin 9 [A 2 + B 2 - C 2 + sin 2 9 (C 2 - A 2 cos 2 - B 2 sin 2 <j>)] }' 
I [2 (C 2 cos 2 9 + B 2 cos 2 sin 2 6 + A 2 sin 2 sin 2 6)] 2 , 



(34) 



A-4cp ab 



= sin 3 9 \ [A 2 B 2 sin 2 9 + C 2 cos 2 9 (A 2 cos 2 <j) + B 2 sin 2 0)] sin 2 (2^) (.BA - AS 



+4 cos 2 9 (B 2 cos 2 <j> + A 2 sin 2 



ACB cos 2 



B ( CA sin 2 



+4 (A 2 - B 2 ) C cos 2 9 cos sin cj> sin(20) ( AB - BA 



AC 
ACB cos 2 



+B (AC - CAsin 2 </>)}/ {2A 2 B 2 [C 2 cos 2 + sin 2 (B 2 cos 2 + A 2 sin 2 </>)] 2 | 



(35) 
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2sin 2 6>|ACcos(6l-0) [-C 2 cos0sin(0 - 0) + vl 2 sin 2 6» sin 0] S 
+B 3 cos sin 2 6» sin(0 - <j>) (cA - AC^j + B cos(6> - 



C 3 Acos0sm( 



-A 3 C sin 2 9 sin < 



} 2 /{{^ 2 [C 2 cos 2 (0-0) + 5 2 sin 2 0] 



and the factor 



-B 2 C 2 sin 2 (6> - <jy)} 2 [C 2 cos 2 + sin 2 (B 2 cos 2 (f> + A 2 sin 2 0)] } , (36) 



/ = e A - 2z = i {A 2 [C 2 cos 2 (0 -<j>) + B 2 sin 2 0] + B 2 C 2 sin 2 (0 - <j>)} . (37) 



V. SIGNIFICANCE OF THE MIXMASTER-[/(l) CORRESPONDANCE 

As was described in Q, the key to the dynamics of the U(l) models is the behavior of the variable ip. In the 
approach to the singularity, one BKL scale factor always dominates. The nature of Mixmaster dynamics is that a 
bounce occurs when one scale factor is w 1 (with the corresponding LSF sa 0). The other scale factors are << 1 with 
the corresponding LSF's large and negative. We shall thus refer to scale factors and functions thereof as "order unity" 
(ps 1) or "exponentially small" (<< 1). During the Kasner epochs (away from the bounces), all the scale factors are 
exponentially small. Assume that during a Kasner epoch A > B > C and define b = B/A, c = C/A, and c = C/B 
so that b, c, and c are < 1. The following discussion can be changed to reflect any other ordering of the scale factors 
by performing the appropriate permutation while keeping in mind that different scale factors will be associated with 



different (0, 4>) dependences. From Eq. (18), we see that tp is dominated by the largest LSF. For our chosen ordering, 



<p = a+ -In [sin 2 9 sin 2 <j> + b 2 (sin 2 8 cos 2 + c 2 cos 2 0)] . (38) 

Since the (ersatz) spatial dependence appears in the logarithm, the spatially homogeneous behavior dominates at any 
particular spatial point. Note also that if b and c are exponentially small, 

Lp ps a + — ln(sin 2 sin 2 (p) w a (39) 

at a generic spatial point. Figure 2 shows a superposition of <p constructed from Eq. Jig) o n the same graph as Fig. 1. 
The scale factors are obtained from a numerical simulation of the Mixmaster ODE's[pl|. Note that ip changes sign 
in two ways. If (p, T > (v v < 0), ip will "bounce" when the scale factor which dominates it bounces. This is the usual 
Mixmaster bounce. However, <p„ also changes sign, after the usual bounce, when the now-increasing scale factor ( B 
in our case) becomes as large as the decreasing dominant scale factor (A in our case). Evaluating V\ from Eq. ( |34| ) 
when A > B > C, we find, for &, c, and c exponentially small, that 

In Vi w 4a + - ln(sin 2 6), (40) 

i.e. bounces off V\ correspond to the usual Mixmaster bounces. This is not surprising since V\ comes from terms in 
the spatial scalar curvature. (It appears in the "kinetic" part of the Hamiltonian constraint because the canonical 
transformation (e a , [3 a ) — * (r, u>) interchanges momentum and configuration variables.) 

If we assume the Kasner time dependence — A — k^A, B — kgB, and C = kcC (recalling that ' = d/dr while the 
Kasner solution is A = A e kAT , for kj\ a constant, etc.), we find in evaluating V2 from Eq. ( |35| ) that the dominant 
scale factor A does not appear so that (for c exponentially small), 

2 (sin + tr cos^ q>) 

implying that V2 is exponentially small unless b = e^~ a is order unity rather than exponentially small. But b = 1 
corresponds to the change in the dominant scale factor from A to B as A decreases while B increases after the standard 



7 



Mixmaster bounce. Thus bounces off V2 in terms of the £7(1) variables correspond to the change in the identity of 
the dominant scale factor as seen in the behavior of ip in Fig. 2. Figure 3 illustrates the interplay among V\_, V2, and 
ip which is very similar to that seen in generic U(l) models (see Figs. 2-5 in |ffl]). 

An effect not mentioned in our previous studies of generic U(l) models appears when we consider the variable z. 
For A > B > C, and b and c exponentially small but c not necessarily << 1, we find from Eq. ( p2] ) that 

sin 2 9\ sin 01 . , 

z « 7 - C + In — 5 ! - (42) 

' S sin 2 9 + ~c 2 cos 2 (0-0) 

— i.e. z decreases monotonically (is dominated by the most negative LSF) unless the era ends and C grows to be as 
large as B so that c w 1. The behavior of z is shown in Fig. 4 for the same simulation as in the other figures. Figure 
5 shows a, £, 7, and z in the region of a bounce which changes the sign of v z . The bounce rules ( [l6|) imply that such 
a sign change is possible if v v > v z /2 + v^/A. Of course, v z < will eventually cause a bounce off V3. 

This rise of z, qualitatively indicating the end of a Mixmaster era, has not been observed in previously published 
simulations of U(l) models Q. However, in Fig. 6a of Q, it is seen that v z changes during bounces off Vi- It is also 
seen that the magnitude of v z decreases at such bounces just as in Fig. 4. Presumably, if the simulations could be 
followed to larger r, z would eventually begin to increase as v z changes sign. Fig. 6a of Q also shows A decreasing 
monotonically with decreasing magnitude of the slope i>a- The bounce rules ( |l6| ) do not permit the sign of v\ to 
change since Eq. (|lj) may be used to show that 6p\ + 2p + p z is always > 0. 

In ifTiyfl, ^ was argued that all terms in H (from Eq. ([l3])) containing spatial derivatives, ex- 
cept for V2, were proportional to / = e A ~ 2z . If A > B > C in a Mixmaster model, / « 
\A A [c 2 cos 2 (6> — <p) + b 2 sin 2 9 + b 2 c 2 sin 2 (6' — <j))] which is always exponentially small since Ami and 6, c << 1. 



Of course, e A = A 4 6esin0y c 2 cos 2 9 + sin 2 9(b 2 cos 2 (f> + sin 2 4>) which is always exponentially small. Thus, we have 
shown that all the observed generic U(l) behavior described in [jjj is characteristic of Mixmaster behavior expressed 
in U(l) variables. 

In order for this to be correct, we also expect the remaining variables u> and x and all the momenta to be of order 
unity at all times so that, generically, they may be regarded as approximately constant in time. This assumption 
underlies the MCP. We find that u> is independent of the scale factors with 

lj = —kA — fcs + sin 2 9(—kc + fcs cos 2 <p + kA sin 2 </>). (43) 

We consider the approximate forms of the remaining variables for the case A > B > C with b, c, and c exponentially 
small. We find 

ikI, p^4kASm9, p z w 2(kc — &b) sin^, 

, s . n 1 (ItB — kc) cos(9 — (h) ,,,, 

p A « -(2k A + k B + k c )sm9, rRJ Px m^ i ZL. (44) 

sm sm sm 



VI. CONCLUSIONS 



Describing a Mixmaster universe as a U(l) symmetric cosmology has yielded significant insight into the observed 
behavior of generic U(l) models. The oscillations in ip observed in the generic models are precisely what would be 
expected from local Mixmaster dynamics. Perhaps of greater interest, a new effect — the signature for the end of an 
era as seen in the rise of z — has been predicted. Observation of this effect in longer running simulations of generic 
U(l) models or in models with an era change at some spatial point early in the simulation would provide strong 
evidence that the observed oscillations are indeed due to local Mixmaster dynamics. Because the rise in z is expected 
to be dramatic, any numerical uncertainties (see the discussion of numerics in Q) should be irrelevant. Simulations 
to look for this effect are in progress. 

One question which remains in our larger program to analyze the approach to the singularity is the dependence 
of the MCP picture and our interpretation of the results on our choice of spacetime slicing and variables. While the 
U(l) description of a Bianchi type IX Mixmaster model uses the same slicing so that no conclusions may be drawn 
on that issue, we clearly have two completely different descriptions of the same Mixmaster model. Both descriptions 
are characterized by intermittent Kasner (or VTD) behavior with bounces off exponential potentials changing one 
Kasner solution into another. However, the three traditional BKL LSF's are replaced in the U(l) description by a 
single variable tp whose oscillations follow the change from one Kasner epoch to the next. The U(l) variable z clearly 
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signals an era change — which is not obvious when following the usual BKL LSF's. In the study of generic collapse, 
we conjecture that, in any convenient variables which allow detection of local VTD behavior, the MCP will predict 
correctly whether or not the model is AVTD. In addition, it should also be possible to identify local Mixmaster 
dynamics through the comparison of the departures from AVTD behavior with the description of a homogeneous 
Mixmaster universe in the same variables. 



APPENDIX 



Using the evolutions equations for the U(l) variables, we find the corresponding momenta to be 

A sin 2 (6>) sin 2 (<p)A + B cos 2 (</>) sin 2 (9) B + C cos 2 (9) C 



p = 4 sin ( 



C 2 cos 2 {9) + sin 2 (6») (B 2 cos 2 (0) + A 2 sin 2 (0)) 



p z = { { -2 sin 9 [{B 2 sin 2 9 cos 2 + C 2 cos 2 9) (A 2 B 2 + A 2 C 2 - B 2 C 2 
-A 2 B 2 cos(20) +(A 2 + B 2 )C 2 cos[2(0 - </>)]) 

-4A 2 B 2 C 2 sin 2 9 sin 2 (9 - 0) sin 2 </>] } j 

+ {2 sin0 [(C 2 cos 2 9 + A 2 sin 2 9 sin 2 0) (-A 2 B 2 + A 2 C 2 - B 2 C 2 
+A 2 B 2 cos(20) + (A 2 + B 2 )C 2 cos[2(0 - </>)]) 

+AA 2 B 2 C 2 sin 2 9 cos 2 (9 - 4>) cos 2 cj>] } | 

+ {8A 2 B 2 C 2 cos 2 9 + (-A 2 B 2 + A 2 C 2 + B 2 C 2 + A 2 B 2 cos(26-) 
+ (A 2 - B 2 )C 2 cos[2(0 - <j>)]) [-A 2 - B 2 

+ (A 2 - B 2 ) cos(20)] } sin 3 0^ j / { [A 2 B 2 + A 2 C 2 + B 2 C 2 

-A 2 B 2 cos(20) + (A 2 - B 2 )C 2 cos[2(0 - </>)] [C 2 cos 2 9 
+ sin 2 9(B 2 cos 2 <j) + A 2 sin 2 0)] } , 

p x = sin 6> {esc 61 [AC cos{6 - <p) (-C 2 cosO sin(0 - <j>) 

+ A 2 sin 2 9 sin <f>) B + B 3 cos sin 2 9 sm{9 - 0) ((J A- AC^j 

+B cos(6> - <j>) (c 3 cos 9 sin(0 - </>) A - A 3 sin 2 9 sin <f> c) } 
[C 2 cos 2 6» + sin 2 6> (B 2 cos 2 + ^ 2 sin 2 </>)]}, 

p A = -sin 9 [2AB 2 Ccos 2 ^sin 2 6lB + ^C (C* 2 cos 2 9 

+ A 2 sin 2 9 sin 2 </>) B + B 3 cos 2 sin 2 (ci + AC 1 ) +S (c 3 cos 2 6>i 

+ 2 A 2 C sin 2 6» sin 2 A + 2 A C 2 cos 2 9 C + A 3 sin 2 9 sin 2 c) 
/{ABC [C 2 cos 2 6> + sin 2 9 (B 2 cos 2 + A 2 sin 2 0)] } . 



(45) 



(46) 



(47) 



(48) 
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FIGURE CAPTIONS 



Figure 1. The LSF's for a portion of a typical Mixmaster trajectory. Note that, prior to w 2, a and 7 oscillate 
while ( decreases monotonically. This "era" ends when £ starts to increase. Now 7 decreases monotonically while a 
and C oscillate. It is clear that another era has ended (but off the scale) near ps 3.5 since 7 is seen to oscillate 
again. 

Figure 2. The U(l) symmetric model variable ip superposed on the Mixmaster trajectory segment of Fig. 1. It is 
clear that <p always tracks the largest LSF and that tp, T changes sign both at the original "bounces" of the LSF's and 
when a different scale factor becomes the largest. 

Figure 3. The U(l) symmetric model variable (p and potentials V\ and V2 for part of the Mixmaster trajectory 
segment of Fig. 1. 

Figure 4. The U(l) symmetric model variable z superposed on the Mixmaster trajectory segment of Fig. 1. Note 
that z decreases monotonically except at the end of an era when it increases. 

Figure 5. Details of the rise of z in Fig. 4. (a) A standard LSF bounce (in () causes z, T to become more negative. 
This is a bounce off V\. (b) The next bounce in <p, off V2 occurs when a — £ and causes z, T to become positive, (c) 
The next bounce (in a) causes z, T to increase, (d) z reaches a maximum when the now increasing smallest LSF 7 
equals the middle one (. 
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